Submitted  to  Applied  Physics  Letters 


Thermal  Conduction  Inhomogeneity  of  Nanocrystalline  Diamond 
Films  by  Dual-Side  Thermoreflectance 


Elah  Bozorg-Grayeli1a),  Mehdi  Asheghi1,  Vincent  Gambin2,  Rajinder 
Sandhu2,  Tatyana  I.  Feygelson3,  Bradford  B.  Pate4,  Karl  Hobart4,  and 
Kenneth  E.  Goodson1 

'Department  of  Mechanical  Engineering,  Stanford  University,  Stanford,  CA,  94305 
2Northrop  Grumman  Corporation,  San  Diego,  CA,  92127 
3Science  Applications  International  Corporation,  McLean,  VA,  22102 
4Naval  Research  Laboratory,  Washington,  DC,  20375 


Abstract 

Thin  diamond  films  of  thickness  near  one  micrometer  can  have  highly  nonuniform  thermal 
conductivities  owing  to  spatially- varying  disorder  associated  with  nucleation  and  grain 
coalescence.  Here  we  extract  the  nonuniformity  for  nanocrystalline  CVD  diamond  films  of 
thickness  0.5,  1.0,  and  5.6  pm  using  picosecond  thermoreflectance  from  both  the  top  and  bottom 
diamond  surfaces,  enabled  by  etching  a  window  in  the  silicon  substrate.  The  thermal 
conductivities  vary  from  less  than  100  W  m  1  K  1  to  more  than  1300  W  m  1  K  1  and  suggest  that 
the  most  defective  material  is  confined  to  within  one  micrometer  of  the  growth  surface. 
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Diamond  is  viewed  as  a  promising  potential  successor  to  Si  and  SiC  as  a  thermal  spreader 
for  high  electron  mobility  transistors  (HEMT)  based  on  GaN1'5.  Since  the  1980’s,  chemical 
vapor  deposition  of  diamond  has  become  well  established6  1 1  and  has  led  to  routine  deposition  of 
polycrystalline  and  nanocrystalline  diamond  thin  films.  A  remaining  challenge  is  to  integrate 
diamond  materials  into  high  power  device  structures  in  order  to  take  advantage  the  high  thermal 
conductivity.  Diamond  thin  film  deposition  requires  a  seed  layer  to  nucleate  film  growth  on  a 
foreign  substrate.  In  this  work,  the  seeds  are  nanodiamonds  (typically  5-10  nm  diamond 
particles)  spread  at  roughly  10  “  seeds/cm"  on  the  silicon  substrate  before  growth.  During  the 
initial  stage  of  deposition,  the  seed  layer  coalesces  into  a  continuous  film.  The  deposited  film  is 
fully  coalesced  at  300nm  thickness.  Beyond  the  initial  coalescence  layer,  there  is  predominantly 
columnar  growth  normal  to  the  silicon-diamond  interface11.  Unfortunately,  if  the  thermal 
resistance  due  to  the  film  coalescence  layer  is  sufficiently  large,  diamond  thin  films  with  a 
coalescence  layer  between  the  bulk  diamond  film  and  the  heat  source  may  not  offer  a  thermal 
advantage  over  Si  or  SiC19  20. 

Theoretically,  isotopically-enriched  single-crystal  diamond  might  eventually  offer 
thermal  conductivities  up  to  5000  W  m  1  K 1  at  room  temperature,  and  experiments  on 
isotopically-enriched  diamond  have  shown  thermal  conductivity  as  high  as  3300  W  m"1  K  1  12  14. 
For  polycrystalline  diamond  films  less  than  several  hundred  micrometers  thick,  the  cross-plane 
thermal  conductivity  remains  near  1000  W  m  1  K  1  at  room  temperature15,16.  Several  factors 
impede  the  mean  free  path  of  phonons  in  these  polycrystalline  diamond  films,  reducing  the 
intrinsic  thermal  conductivity.  These  include  impurity  scattering,  grain  boundary  scattering,  and 

i  c  i  o 

interface  scattering  '  .At  the  initial  stage  of  diamond  thin  film  deposition,  these  scattering  sites 
are  abundant,  significantly  reducing  thermal  conductivity  .  This  coalescence  region  is  buried 
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underneath  several  jam  of  columnar  diamond  grains,  making  the  thermal  properties  difficult  to 
access.  Since  the  coalescence  region  strongly  affects  the  diamond  thermal  performance,  it  is 
critical  to  develop  a  technique  which  can  extract  the  thermal  properties  of  both  the  coalescence 
and  columnar  diamond  region. 

Here  we  perform  picosecond  time-domain  thermoreflectance  (TDTR)  measurements  on 
suspended  nanocrystalline  diamond  films  from  0.5  pm  to  5.6  pm  in  thickness.  By  probing  both 
sides  of  the  suspended  film,  we  capture  the  cross-plane  thermal  conductivity  of  both  the 
coalescence  region  ( kc )  and  the  high-quality  columnar  grain  (Jchq)  regions.  Using  these  results 
with  a  two-layer  heat  diffusion  model  of  the  diamond  film,  we  estimate  the  thickness  of  the  low 
conductivity  coalescence  layer  and  compare  the  thermal  resistances  of  the  coalescence  versus 
high-quality  regions.  This  analysis  demonstrates  that  the  coalescence  region  can  be  a  substantial 
contributor  to  thermal  resistance  in  a  diamond  thin  film. 

Picosecond  TDTR  uses  the  temperature-dependent  optical  properties  of  a  thin  metal 

21  22 

transducer  as  a  thermometer  for  a  heated  thin  film  stack  ’  .  Our  TDTR  system  uses  a  passively- 
modelocked  Nd:YV04  laser  with  9.2  ps  pulsewidth  and  82  MHz  repetition  rate"  ’  .  We  then  use 
a  radially- symmetric  model  of  heat  diffusion  through  a  multilayer  stack  to  fit  for  the  transducer- 
film  thermal  boundary  resistance  (TBR),  the  intrinsic  thermal  conductivity  of  the  film,  and  the 
film- substrate  TBR26'27. 

The  heat  capacitance  of  the  diamond  coalescence  region  is  small  compared  to  the  total 
heat  capacitance  of  the  diamond  film,  and  the  coalescence  region  is  several  pm  from  the 
transducer  film.  As  a  result,  the  additional  coalescence  layer  thermal  resistance  is  impossible  to 
extract  separately  from  the  diamond-silicon  thermal  boundary  resistance"  .  In  order  to  separate 
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the  two  properties,  we  require  direct  contact  between  the  transducer  layer  and  the  coalescence 
layer.  This  involves  removal  of  the  silicon  substrate  to  allow  deposition  of  a  metal  transducer 
directly  on  the  seeded  surface  of  the  coalescence  region.  To  accomplish  this,  we  deposited  three 
nanocrystalline  diamond  films  on  silicon  substrates  of  thicknesses  0.5  pm,  1.0  pm,  and  5.6  pm. 
The  silicon  substrates  were  selectively  etched  to  produce  a  suspended  diamond  thin  film  of 
roughly  5  pm  in  diameter.  A  50  nm  A1  transducer  layer  was  deposited  on  each  side  of  the 
suspended  diamond  film  (fig.  1). 
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FIG.  1.  Cross-sectional  diagram  of  the  suspended  diamond  fdm  sample  showing  the  Al  transducer  interface  with  the 
seeding/coalescence  layer  and  with  the  high-quality  (HQ)  nanocrystalline  diamond  surface.  Note  that  silicon  has 
been  etched  to  give  access  to  the  seeded  surface  of  the  diamond  fdm  where  fdm  coalescence  occurs.  The  Al 
transducer  layer  is  50  nm  thick,  and  the  etched  hole  in  silicon  is  ~  5  mm  in  diameter.  The  red  and  green  areas  (color 
image  available  online)  indicate  the  probe  and  pump  beams,  respectively. 


When  the  characteristic  thermal  decay  time  of  the  film  of  interest  is  less  than  the 


measurement  time  scale,  and  the  film  is  effectively  insulated  during  the  measurement  time  scale, 


the  temperature  decay  of  the  transducer  film  depends  primarily  on  the  heat  capacity  and  total 


thermal  resistance  of  the  film  of  interest.  For  the  case  of  the  0.5  pm  film,  we  obtain  only  an 

2  1 

effective  thermal  resistance  for  the  diamond  film  of  20  ±  2.0  m  K  GW’  ,  which  can  be  expressed 


as: 
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where  Rey  is  the  effective  resistance  of  the  diamond  film,  RAi-diam  is  the  aluminum-diamond  TBR, 
and  ddiam  is  the  thickness  of  the  diamond  film.  We  convert  this  to  an  effective  thermal 
conductivity  of  the  diamond  layer  by  dividing  the  film  thickness  by  the  resistance.  This  gives  a 
value  of  keff,c  of  approximately  25  W  m  1  K1.  Using  the  previously  reported  diamond  heat 
capacity  is  ~  2.15x10  J  m"  K  ,  this  translates  to  a  total  thermal  diffusivity  of  ~  1.16x10"  m“ 
s’1.  The  characteristic  thermal  decay  time  of  the  diamond,  determined  via: 


T-diam 


d2 

** diam 


a, 


eff 


(2) 


where  aey  is  the  effective  thermal  diffusivity,  is  ~  22  ns.  The  characteristic  timescale  of  the 
heating  event,  however,  is  given  by  the  inverse  of  the  5  MHz  pump  modulation  frequency.  This 
translates  to  a  heating  timescale  of  200  ns,  significantly  greater  than  Xdiam-  Therefore,  although 
measurements  on  this  sample  are  insensitive  to  the  difference  between  intrinsic  and  thermal 
boundary  resistance,  we  can  uniquely  extract  the  diamond  heat  capacity.  The  fitted  heat  capacity, 
1.98xl06  J  nr  K'1,  is  similar  to  previous  literature  results28. 

In  order  to  separate  out  kc  and  RAi-diam  from  the  effective  thermal  resistance,  we  use  the 
1.0  pm  thick  diamond  film.  While  for  this  film  (~56  ns)  is  still  less  than  the  timescale  of  the 
measurement,  we  have  sufficient  sensitivity  to  separate  the  transducer-film  TBR  from  the 
intrinsic  film  thermal  conductivity.  Modeling  the  system  as  50  nm  A1  on  1.0  pm  diamond,  we 
extract  a  TBR  of  13.5  ±  1.0  nr  K  GW"  between  the  A1  transducer  and  the  diamond  coalescence 
layer.  Further,  we  find  a  cross-plane  thermal  conductivity  of  80  ±  10  W  m  1  K  1  for  the  diamond 
coalescence  layer.  By  applying  these  results  to  equation  (1),  and  setting  ddum  to  0.5  pm,  we 
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obtain  a  total  thermal  resistance  that  agrees  with  the  value  measured  from  the  0.5  pm  diamond 
sample.  Comparing  the  results  with  the  model  of  local  diamond  thermal  conductivity  developed 
by  Touzelbaev  et  al  indicates  an  average  coalescence  layer  grain  size  of  100-200  nm  . 

We  determine  the  thickness  of  the  thermally  resistive  coalescence  layer  using  the  5.6  pm 
diamond  film.  Picosecond  TDTR  measurements  on  the  top  side  of  the  diamond  reveal  an  Al- 
diamond  thermal  boundary  resistance  of  10.7  ±  1.0  m  K  GW"  ,  smaller  than  for  the  Al- 
coalescence  layer  interface.  Since  the  thermal  boundary  resistance  is  highly  dependent  on  surface 
cleanliness,  this  difference  may  be  due  to  impurities  left  behind  by  the  etching  process  on  the 
bottom  diamond  film  surface.  Assuming  the  heat  capacity  measured  from  the  0.5  pm  sample,  and 
using  a  one-layer  model  of  50  nm  Al  on  5.6  pm  diamond,  we  find  a  cross-plane  thermal 
conductivity  of  1350  ±  200  W  m’1  K  1  for  the  high-quality  diamond  (fig.  2).  It  is  worth  noting 
that  the  result  of  the  numerical  fit  for  this  sample  remains  the  same  even  assuming  a  diamond 
thickness  as  low  as  2.5  pm.  This  implies  the  top  side  measurement  may  be  insensitive  to  the 
coalescence  layer  properties. 

To  confirm  this  hypothesis,  we  performed  TDTR  on  the  bottom  of  the  same  sample. 
Figure  3  demonstrates  the  difference  in  the  thermoreflectance  curves  for  the  two  measurements, 
showing  a  significantly  slower  thermal  decay  through  the  coalescence-layer.  Using  the 
coalescence  layer  diamond  data  from  the  5.6pm  thick  film,  we  created  a  two-layer  model  of  heat 
conduction  through  the  suspended  film.  This  model  assumes  50  nm  of  Al  on  a  diamond 
coalescence  layer  with  unknown  thickness,  dc,  on  a  high-quality  diamond  layer  of  thickness  dHQ 
=  5.6  pm -dc-  Assuming  kc  =  80  W  m"1  K  1  and  kHQ  =  1350  W  m"1  K'1,  we  fit  for  the  coalescence 
layer  thickness  and  find  dc  =  0.76  ±0.1  pm.  The  thickness  of  the  high  quality  diamond  layer, 
therefore,  is  -4.84  pm,  thick  enough  to  render  the  top  side  measurements  insensitive  to  the 
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thermal  properties  of  the  coalescence  layer  diamond.  This  validates  the  one-layer  assumption 
used  for  the  top  side  measurement  of  the  5.6  pm  sample. 


FIG.  2.  Picosecond  TDTR  data  (solid  line)  and  numerical  fit  (dotted  line)  for  the  top  side  measurement  of  the  5.6 
pm  diamond  film.  The  fit  corresponds  to  an  Al-diamond  thermal  boundary  resistance  of  10.7  ±  1.0  m2  K  GW'1  and  a 
diamond  thermal  conductivity  of  1350  +  200  W  m'1  K1. 


FIG.  3.  Comparison  of  picosecond  TDTR  curves  from  the  coalescence  and  high  quality  layers  of  the  5.6  pm 
diamond  sample.  The  larger,  higher  quality  diamond  grains  on  the  top  surface  result  in  a  higher  thermal  conductivity 
than  for  the  bottom.  This  is  evidenced  by  the  faster  thermal  decay  shown  for  the  high  quality  layer. 

The  thickness  prediction  for  the  coalescence  layer  is  more  subtle  than  the  two-layer 
model  implies.  Although  the  diamond  film  consists  of  columnar  grains  on  a  thin  coalescence 
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layer,  there  is  no  sharp  transition  between  the  two.  Rather,  the  average  grain  size  in  the  film 
increases  smoothly  as  a  function  of  distance  from  the  substrate  surface.  What  the  two-layer 
model  determines  is  an  estimated  thickness  based  on  the  thermal  conductivity  measured  in  the 
thinner  samples.  If  we  assume  a  higher  thermal  conductivity,  the  fitted  coalescence  layer 
thickness  increases  proportionally.  A  more  robust  comparison  would  involve  the  thermal 
resistances  of  the  two  layers.  Although  the  fitted  thickness  of  the  coalescence  layer  depends 
proportionally  on  the  assumed  thermal  conductivity,  the  fitted  thermal  resistance  of  the  layer 
remains  the  same  regardless  of  assumed  thermal  conductivity. 

The  thermal  properties  of  the  diamond  film  coalescence  layer  have  a  drastic  effect  on  the 
total  thermal  resistance  of  the  film.  Using  the  coalescence  layer  thickness  and  thermal 
conductivity  we  estimated  using  the  two-layer  model,  we  find,  for  a  5.6  pm  diamond  film,  that 
the  coalescence  layer  thermal  resistance  ( Rc  —  dc/kc)  is  ~  9.5  ±  1.4  m2  K  GW"1.  This  is 
significantly  larger  than  the  thermal  resistance  of  the  rest  of  the  diamond  film  (~  3.6  ±  0.6  m“  K 
GW"1  assuming  a  0.76  pm  coalescence  layer).  In  fact,  the  coalescence  layer  has  a  thermal 
resistance  equivalent  to  ~  12.8  pm  of  high  quality  diamond. 

This  work  extracts  the  thermal  conductivity  nonuniformity  in  nanocrystalline  diamond 
films  of  thickness  near  5  micrometers  using  dual-sided  picosecond  TDTR  measurements.  This 
technique  extracted  the  thermal  conductivity  of  the  coalescence  layer,  the  high  quality  diamond 
layer,  and  the  heat  capacity  of  the  suspended  films.  From  these  results,  we  estimate  the 
coalescence  layer  thickness  and  calculate  the  thermal  resistance  contributions  of  both  regions  of 
the  diamond  film.  The  additional  coalescence  layer  thermal  resistance  is  significant.  Further 
research  into  nanocrystalline  diamond  nucleation  and  growth  may  improve  the  coalescence  layer 
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thermal  properties,  reducing  the  total  thermal  resistance  and  enhancing  the  thermal  performance 
of  integrated  nanocrystalline  diamond  thermal  spreaders  for  high-power  transistor  applications. 
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